Dipole strength in Sm studied via (7,n), (7,p) and (7, a) reactions 
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Photoactivation measurements on ^^''Sm have been performed with bremsstrahlung endpoint ener- 
gies from 10.0 to 15.5 MeV at the bremsstrahlung facility of the superconducting electron accelerator 
ELBE of Forschungszentrum Dresden-Rossendorf. The measured activation yield for the ^**Sm(7,n) 
reaction is compared with the calculated yield using cross sections from previous photoneutron ex- 
periments. The activation yields measured for all disintegration channels ^^*Sm(7, n), (7,p) and 
(7, q) are compared to the yield calculated by using Hauser-Feshbach statistical models. A new 
parametrization of the photon strength function is presented and the yield simulated by using the 
modified photon strength parameters are compared to the experimental data. 

PACS numbers: 25.20.-x, 25.20.Dc, 26.30.-k 
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I. INTRODUCTION 

The nuclei heavier than iron [Z > 26) are synthesized 
mainly by neutron-capture reactions- the astrophysical 
r- and s-processes. There are about 35 neutron deficient 
stable isotopes between Se and Hg that are shielded from 
the rapid neutron capture by stable isobars. These nu- 
clei, classically referred to as the p-nuclei, are produced 
via chains of photodisintegrations such as (7,^) (7,73) 
and (7, a) on r- or s-seed nuclei [l|, [2]- The best possi- 
ble sites that have been proposed for the production of 
p-nuclei are the 0/Ne-rich layers of type II supernova 
explosions [3]. 

For the calculation of p-process abundances, a compre- 
hensive network involving thousands of reaction rates is 
necessary. Despite the efforts in recent years, the experi- 
mental information on the reaction rates involved in the 
p-process nuclear flows is very scarce. The p-process re- 
action rates presently used in the astrophysical network 
calculations are based on the cross sections obtained from 
Hauser-Feshbach statistical model calculations. 

For the network calculations involving the p-nuclei, the 
precise knowledge of the photon strength especially in the 
close-threshold regions is very important. In contrast to 
the rather detailed knowledge [4| of the photon strength 
in the isovector giant dipole resonance (GDR) region well 
above the particle separation energies, the corresponding 
information for the excitation region near threshold is 
still surprisingly uncertain due to the much smaller cross 
section. Accurate information on weaker channels with 
protons or alpha particles in the exit channels is also 



'Present address: Nuclear Engineering Division, Argonne National 
Laboratory, Argonne, Illinois 60439, USA. 

t Present address: Istituto Nazionale di Fisica Nucleare, Sezione di 
Padova, 35131 Padova PD, Italy. 

-^Present address: Department of Physics, Duke University and Tri- 
angle Universities Nuclear Laboratory, Durham, NC 27708, USA. 



necessary. 



The experimental information available for the p- 
nucleosynthesis so far is either from the reactions that 
involve neutrons or charged particles. The work de- 
scribed here presents experimental data for the photo- 
disintegration of one nucleus into four different chan- 
nels as observed in a wide range of bremsstrahlung end- 
point energies, i.e, the ^*^Sm(7, n) to ^^^^Sm and ^''•^'"Sm, 
the "'^*^Sm(7,p) and ^^^Sm(7,a) reactions. The nuclide 
^*^Sm has been discussed in the frame of the p-process 
chronometer ^"^^Sm ^, 61]. There were several efforts to 
determine the ^"'^Sm/^''^Sm production ratio experimen- 
tally which varies due to uncertainties in different inputs 
entering into the calculation, see Ref. [T]. The nuclear- 
physics-related uncertainties in the p-process modelpre- 
dictions were investigated recently by Rapp et al. [8|. A 
list of the most critical p-process reaction rates (with 
their respective inverse reactions), which will influence 
the final p-abundances, have been given (see Table 2 and 
Table 3, Ref. Q). In general, a systematic investigation 
of the (7, a) reactions in the mass region A>140 is found 
to be very important. Therefore, experimental informa- 
tion on the photodisintegration rates in ^*^Sm helps to 
reduce the inherent nuclear model uncertainties. 



The photoactivation experiments presented here cover 
a substantial range of excitation. The experimental setup 
is described in SecHIland Sec. lIIII explains the method of 
photoactivation yield determination. The experimental 
results of the ^^''Sm(7, n), ^^''Sm(7,p) and ^^^Sm(7, a) 
reactions are discussed in detail under Sec. IIVI The pho- 
toactivation experiments are performed such that an ab- 
solute yield can be extracted which allows the derivation 
of an absolute strength. This strength function is used 
in statistical calculations of Hauser-Feshbach type as an 
input parameter. This point is worked out in detail in 
Sec. lIVBl of this paper. 



II. EXPERIMENTAL SETUP 

The photodisintegration experiments discussed here 
were performed at the bremsstrahlung facihty of the su- 
perconducting electron accelerator ELBE (Electron Lin- 
ear accelerator of high Brilliance and low Emittance) 
of Forschungszentrum Dresden-Rossendorf. At ELBE, 
bremsstrahlung endpoint energies up to 20 MeV and av- 
erage currents up to 1 mA are available which is ap- 
propriate for probing photon-induced reactions. The 
bremsstrahlung facility has been extensively used for 
photon scattering as well as photoactivation studies [3- 
[l5| . The experimental setup at ELBE has been described 
in detail elsewhere [la, [131 . 

A schematic layout of the photoactivation facility is 
given in Fig. [T] The primary electron beam is focused 
onto a thin radiator foil which produces bremsstrahlung 
via deceleration of electrons. The radiator is made of 
niobium with areal densities varying between 1.7 mg/cm^ 
and 10 mg/cm^ corresponding to 1.6 x 10"'' and 1 x 10~^ 
radiation lengths. Behind the radiator, the electrons 
are deflected by a dipole magnet and dumped into 
a graphite cylinder with a conical recess (see Fig. [l] 
photoactivation site). The length of the cylinder is 600 
mm and diameter is 200 mm. The photoactivation tar- 
get (Sm) is irradiated here together with the activation 
standard target "'^^Au. 

The bremsstrahlung beam goes straight ahead through 
the collimator to the photon scattering site. The photon- 
scattering site is separated from the photoactivation site 
by a 1.6 m thick heavy-concrete wall. The collimator 
made of high-purity aluminum is placed 1 m behind the 
radiator. The collimator shapes a beam with defined di- 
ameter from the spatial distribution of photons. For the 
energy range of interest, the flux at the photon scatter- 
ing site is about 10^ cm~^ s~^ MeV~^. The endpoint 
energy of the bremsstrahlung distribution is determined 
by measuring protons from the photodisintegration of the 
deuteron (see Fig. [I] deuteron breakup target) with sili- 
con detectors. From the maximum energy of the emitted 
protons, the maximum energy of the incident photons 
was deduced. For a detailed description of the beam en- 
ergy determination procedure, see Sec. Ill C of Ref. [l5| . 

At the photon scattering site, a "'"'B sample is irradi- 
ated together with the activation standard target ^^^Au. 
The activation standard reaction ^^^Au(7, n) has been 
compared to photoneutron studies using monochromatic 
photons from positron annihilation in flight technique 
(see [l3|, and the references therein). The photon flux 
at this site is experimentally determined by means of the 
known integrated cross sections of the states in ^^B de- 
populating via 7 rays. The thin target bremsstrahlung 
and photon flux determination procedure has been dis- 
cussed in detail in Ref. 15]. 

At the photoactivation site, the available photon flux 
is about 50 to 100 times higher compared to the photon 
scattering site. At this high-flux area, it is technically 
not possible to determine the photon flux online. The 
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FIG. 1: (Color online) Photoactivation setup at ELBE ac- 
celerator [1^. The electron beam is deflected from the main 
beam line and creates bremsstrahlung in the radiator. There 
are two target sites: at the photoactivation site, Sm target is 
irradiated together with ^®^Au which is used as a reference. 
At the photon scattering site, the scattered photons from ^^B 
are observed using HPGe detectors. Another ^®^Au reference 
target is sandwiched with ^^B for flux calibration purposes. 



bremsstrahlung distribution by the graphite is described 
by MCNP4C2 pLg simulations which are based on the 
bremsstrahlung cross sections by Seltzer and Berger [l9| . 

For the activation measurements, the Sm/Au samples 
used were of natural isotopic composition with masses 
^-^200 mg (Au) and ~2 g (Sm). The Sm powder was 
pressed into samples of cylindrical shape and the masses 
were determined prior to the irradiation. After irradi- 
ation, the decay of the daughter nuclei resulting from 
photoactivation is studied with HPGe detectors of rel- 
ative efficiency 90% or 60% situated in a lead-shielded 
low-background environment. Between the sample and 
the endcap of the detector, a cadmium absorber with 
a thickness of 1.535 mm is inserted to block the X-ray 
summing with the 7-rays. To maximize the absolute effi- 
ciency of the detector, the source/target is put on top of 
the Cd absorber which is situated directly on the HPGe 
capsule. 

In case the resulting radioactive nuclei are short-lived, 
a rapid transport system {rabbit system) was used for ac- 
tivation experiments. The rabbit system uses compressed 
air to transport the samples between the photoactivation 
site and the measuring site [20|. The samples to be ir- 
radiated are enclosed in polyethylene cassettes and are 
shot to the high photon flux area (photoactivation site. 
Fig. [1]). After irradiation, the samples are transported 
within about 15 s to the detector for decay measure- 
ments. For the evaluation of fast decaying samples, a 
list mode data acquisition system has been used which 
enabled us to account for the rapidly changing dead time 
fractions during the course of the measurement. 

The measurements demanding low background 
were performed at the underground laboratory 
'Felsenkeller' [21[ in Dresden where 98% of the cosmic 
niuons are shielded by a 47 m thick rock layer. At 
Felsenkeller, the 7-decay was measured with a HPGe 
detector of 30% relative efficiency. 



III. DATA ANALYSIS 

The data analysis method is explained in the following 
sections. The discussion is divided into two parts: 

(A) the method of photoactivation and the photoactiva- 
tion yield determination; and 

(B) the photodisintegration reactions {'y,n), (7,p) and 
(7, a) in "'^^^Sm. 

The decay properties of the radioactive nuclei resulting 
from photoactivation are given in Table ID 



A. The photoactivation yield 

The photoactivation method essentially consists of two 
steps 

(i) irradiation of the target nuclei ("'^'^^Sm ) 

(ii) measurement of decay of the daughter nuclei result- 
ing from photoactivation 

The irradiation is performed at the photoactivation site 
shown in Fig.[T] The number of activated nuclei A^actl-E-o) 
produced is proportional to the integral of the absolute 
photon fluence (time integrated flux) ^■y{E, Eq) times the 
photodisintegration cross section (T(^_x)(^)- The integral 
runs from the reaction threshold energy Ethr up to the 
endpoint energy Eq of the bremsstrahlung spectrum. The 
symbol x — n,p,a denotes the emitted particle. The 
number of target atoms in the sample is denoted by iVtar- 



iVact(So) = A^tar ' / ^(^.x)(S) ' "^-^iE, Eq) dE (1) 

After irradiation, the 7-rays following the /3-decays of the 
radioactive sample are measured using HPGe detectors. 
The number of activated nuclei N^ctiEo) is determined 
using the relation: 



A^act(-Bo) 



Nj{E^,Eo) ■ Kco 
e{E^)-p{E^) 



(2) 



N^{E^,Eo), e{E^) and p{E^) denote the dead-time and 
pile-up corrected full-energy peak intensity of the ob- 
served transition, the absolute efhciency of the detector 
at the energy E^ and the emission probability of the pho- 
ton with energy Ej respectively. 

The factor Kcorr in Eq. ([2]) is given by 



exp(^ 



l-exp(-^) l-exp(-^ 



(3) 



This expression determines the number of radioactive nu- 
clei from their decays measured during the time tmoas- It 
also takes into account decay losses during irradiation 
(tjir) and in between end of irradiation and beginning of 



measurement (tioss)- t denotes the mean life time of the 
radioactive nucleus formed during photoactivation. 

For the (7,2;) reaction (x = n,p or a), the activation 
yield is denoted by Yact and is expressed as the ratio of 
the number of activated nuclei to the number of target 
atoms in the sample. For example, for the ^^^Sm(7,n) 
reaction. 
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(4) 



Using Eq. ([T]) , the activation yield can be calculated from 
f (7,n) (E) data with the known bremsstrahlung spectrum. 
In this way measured activation yields can be compared 
with the experimental or theoretical cross section data. 

To compare the activation yield measured at different 
endpoint energies with the calculated yield using cross 
sections from theory/previous experiments, the experi- 
mental data need to be normalized to the photon flu- 
ence at the irradiation site (the photoactivation site, see 
Fig.[T]). As already mentioned, it is not possible to mea- 
sure the photon fluence at this site directly. The flu- 
ence at a flxed energy E^ is given by the ratio of the 
measured "'^^^Au(7, n) activation yield and the calculated 
activation yield using the known crij^n) from ^^''Au and 
a simulated thick target bremsstrahlung spectrum using 
the code MCNP. For the experiments discussed in this pa- 
per, the normalization has been done to the fluence at E^ 
= 7.288 MeV which is well below the bremsstrahlung end- 
point energies for the photodisintegration experiments 
under discussion. 

For the precise and accurate analysis of the decay spec- 
tra, the HPGe detectors were calibrated with a certified 
set of calibration sources in the energy range from 0.12 
to 1.9 MeV. The efficiency curve for the detector with 
90% relative efficiency is shown in Fig. [5J For multi- 
gamma emitting nuclides, coincidence summing correc- 
tions have been applied. The fit to the data points is a 
curve based on a GEANT3 [2^ simulation with realistic 
detector geometry. The detector dimensions have been 
determined independently by performing a 200 keV X- 
ray scan. The simulated and measured efficiencies agree 
within 1.01±0.02 for the considered energy range. In the 
estimation of errors, both statistical and systematic un- 
certainties are taken into account. The systematic uncer- 
tainties in the activity as given in the source certificates 
amount to 0.6-1.5% in the energy range from 0.12 to 1.9 
MeV. The statistical uncertainty mainly originates from 
the peak fit for the calculation of peak areas and the 
dead-time and pile-up corrections. The statistical uncer- 
tainty contribution is in the order of 0.2-0.5%. However, 
for the photoactivation experiment described here, the 
targets used were not point-like and the appropriate cor- 
rection for volume-effects has been applied. This will be 
described in the next section. 



TABLE I: Decay properties of the daughter nuclei stemming from the respective photodisintegration reactions in ^^''Sm. The 
values are adopted from the Evaluated Nuclear Structure Data Files of the National Nuclear Data Center online service [23 • 



reaction 



5'.° (MeV) 



decay 



^1/2 



E^" (keV) 



*Sm(7,5 



10.5200(24) 



'Sm(EC+/3+)"^Pm 



^Sm IT decay 



8.75(8) min 


1056.58 


0.019(2) 




1173.18 


0.004161(544) 




1514.98 


0.006593(75) 




1403.06 


0.003496(414) 


66(2) s 


754.4 


0.8988(6) 



*Sm(7,a) 



-0.145 



i40Nd(EC)"°Pr(no 7) 3.37(2) d 

i«Pr(/3+)i«Ce 



1596.11 



0.0050(4) 



'Sm(7,p) 



6.295(3) 



^Pm(EC)'''^Nd 



265(7) d 



741.98 



0.385(24) 



"separation energy for the proton or neutron emission; Q-valuc 
for alpha emission 

''half-life of the corresponding decay with absolute uncertainty 
given in parentheses 

'^7-energy of the analyzed transitions 

''photon emission probability per decay with absolute uncertainty 
given in parentheses 
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FIG. 2: Full-energy peak efficiency curve for a HPGe detector 
with 90% efficiency. The upper part is the fit curve based on 
GEANT3 simulations normalized to the experimental data 
and the lower part shows the relative residuals. 
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FIG. 3: (Color online) "*Sm(7,7') (green),"''Sm(7, n) (blue 
to the ground state, black to the isomer), ^''■'Sm(7,p) (pur- 
ple) and ^*'*Sm(7,a) (red) cross sections from model calcu- 
lations. TALYS and NON-SMOKER predictions are shown 
with dashed and solid lines respectively. The TALYS calcu- 
lations were performed with the default input parameter set 
in the version TALYS-1.0. 



B. Photodisintegration of ^**Sni 



In this experiment, the photodisintegration reactions 
(7, n), (7,p) and (7, a) on ^'*^Sni were studied. The cross 
section predictions for the (7,^), (7,p) and (7, a) reac- 
tions by the statistical models TALYS ^ and NON- 
SMOKER [H are given in Fig. O The i44Sm(7, n) reac- 
tion produces ^^^Sm or "Smg^^ rpj^^ NON-SMOKER 



model doesn't provide the cross section data for the 
ground and isomeric states in ^^^Sm separately. Above 
11 MeV, the (7,p) and (7,0;) cross sections are orders 
of magnitude lower than the (7,^1) cross sections due to 
the suppressive effect of the Coulomb barrier. The re- 
action thresholds and decay properties of the respective 
photodisintegrations in ^'^'^Sm are given in Table HI 

For studying the ^'**Sm(7, n) reaction, the irradiation 



was performed for various endpoint energies starting with 
Eq = 11.00 MeV which is approximately 480 keV above 
the neutron separation energy. The daughter nucleus 
produced is ^^'^Sm, partiahy in its isomeric state ^^■^Sm™. 
The ^■*'*Sm(7, n)^^'^Sm™ reaction is identified with the 
isomeric transition at 754 keV and the ^^'*Sm(7, n)^'*'^SmS 
with transitions above 1000 keV (see Table|l]). Since both 
radionuchdes are short-lived, the irradiation was carried 
out using the rabbit system. The activations for this 
experiment lasted between 15-45 minutes. The sample 
spectrum of an irradiated Sm target and the half-life mea- 
surement using the exponential decay of the ground state 
of ^^-^Sm have been presented in Fig. 2 of Ref. |20| . 

For the activation using the rabbit system, the targets 
used were not point-like. For the irradiation of samarium, 
we used fine Sm203 powder filled in a cylinder (length 16 
mm, diameter 10 mm) of about 2-3 g mass. The ^^^Au 
activation standard targets were made of Au-metal foils 
of about 200-250 mg mass rolled into a cylindrical shape 
(length 36 mm, diameter 5 mm). Therefore, the full- 
energy peak efficiencies measured with point sources have 
to be corrected for source-extension and self-absorption 
effects. To account for the volume source effect, detailed 
simulations using GEANT3 were performed. The abso- 
lute full-energy peak efficiency of the target is given by 



^target [-^'f ) — ^point v-^7 



(E^) 



-point y-^l 



(E^) 



(5) 



£point denotes the efficiencies measured with point-like 
sources at a distance corresponding to the distance of 
the center of the volume source to the HPGe crystal. 
The simulated full-energy peak efficiencies for volume 
and point-like sources are denoted by ^target ^'^^ ep™nt 
respectively. The ratio is the correction to be applied for 
the volume source effect. For example, for the transition 
at 754 keV resulting from the ^^■^™Sm decay, the correc- 
tion factor f^^Hii^) amounts to about 0.97±0.03. For 

the transition at 356 keV following the decay of ^^^Au, 
the correction factor was about 0.93±0.02. 

The ^*^Sm(7,p) and ^'^^Sm(7, a) activation yields were 
measured for the first time at astrophysically relevant en- 
ergies. In the context of this work, these are categorized 
as manual measurements, indicating that the rabbit sys- 
tem was not necessary for the transport of the rather 
long-lived isotopes to the counting setup. 

For the energy range considered in these experiments, 
the ^'*^Sm(7,p) and (7, a) cross sections are rather very 
small compared to the (7,71) cross sections (see Fig. [3|). 
This results in a decay spectrum with weak counting 
statistics. Therefore, the decay measurements were per- 
formed under optimized background conditions in the 
'Felsenkeller' underground laboratory. 

In the photoactivation of ^^'*Sm, the (7,^) and (7,^) 
reaction yields cannot be distinguished at high energies 
above the neutron separation energy (S'„ = 10.52 MeV) 
since the (7, n) daughter nucleus decays very quickly to 
the (7,7?) daughter nucleus. Therefore, to measure the 
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FIG. 4: The spectrum of an irradiated 8111203 target in which 
the decays following the ^*^Sin(7,p) and ^^^Sni(7,Qf) disinte- 
grations are observed. The strongest 7-transitions following 
the /3-decays of ^^'^Pm and ^^''Pr are marked. The target was 
placed on top of the 30% HPGe detector at the underground 
laboratory. The other peaks are from the contamination re- 
actions like ^^''Sm(7, n) etc. The target was irradiated at a 
bremsstrahlung endpoint energy of 15.0 MeV and the spec- 
trum shown here was recorded for a period of 1 day. 



pure (7,p) reactions, we used bremsstrahlung energies 
sufficiently below 10.52 MeV. 

A sample spectrum of an irradiated samarium target is 
given in Fig. S) The Sm203 target was placed on top of 
the 30% HPGe detector. In this measurement, the tar- 
get was irradiated at a bremsstrahlung endpoint energy 
of 15.0 MeV and was brought to the underground labora- 
tory on the next day of irradiation. The spectrum shown 
here was recorded for a period of 1 day. The ^'*^Sm(7,p) 
reaction is identified by the transition at 742 keV from 
^^■^Pm decay (EC) with a half-life of 265 days, see figure. 
The ^^^Sm(7,Q!) reaction was identified by the transition 
at 1596 keV following the decay of ^''Opr which is the 
daughter of the (7, a) product ^'^''Nd. In the spectrum 
shown in Fig. SI the 1596 keV transition is also marked. 
The other peaks in the spectrum stem from the contam- 
inating reactions in the other samarium isotopes (e.g., 
i54Sm(7,n)). 

For the manual measurements, we used Sm2 03 pow- 
der filled into discs of radius 9 mm and thickness 5 mm. 
To correct the full-energy peak efficiencies for the volume 
source effect, we used the Monte-Carlo efficiency trans- 
fer code EFFTRAN [26.]. The point source efficiency 
has been measured with the point source shifted to a 
distance which corresponds to the center of the volume 
source. The point source efficiency is reduced by 7% 
when the point source is shifted and put in the middle of 
the volume target. The volume source effect is only 4%, 
comparing the point source at the center of the Sm203 
sample and a volume target Sm203 placed on the top of 
the cadmium absorber. The ^^^Au activation standard 



targets were very thin discs (masS'^lOO mg, thickness-0.2 
mm, radius-9 mm) and the volume source effect was not 
significant. 

The activation yield relative to the photon fluence are 
presented in the next sections in detail. 



IV. RESULTS AND DISCUSSION 

A. Activation yield for ^''^Sm(7,ri) reaction: 
Comparison to previous experiments 

The activation yield for the ^ Sm(7, n) reaction was 
determined using the method discussed in Sec. lIII Al The 
experimental activation yield for the ^''^Sm(7,n)^^'^Sm 
reaction (Yg) relative to the photon fluence is given in 
Fig. El 

The activation yield is compared to the yield calculated 
using cross sections from previous experiments. The pho- 
toneutron cross sections of the Sm isotopes have been 
measured by Carlos et al. ^25] using positron annihila- 
tion in flight beams at Saclay. In particular, Carlos et al. 
present a study of the transition from spherical to de- 
formed shape for isotopes in the samarium region. The 
partial photoncutron cross sections - [g(7,n)+a(7,np)] 
and (7(7, 2n) of ^'^''Sm are given in Fig. 2 of J27J . 
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FIG. 6: The isomeric yield ratio r] = y " 

ELBE (diamonds) for the ^^*Sni(7,n) reaction is compared 
to the isomeric yield ratio given by Mazur et al. (open dia- 
monds) . 



data obtained with natural Sm at the same labo- 
ratory [28| at the same photon energies. For the 
energy range 11-16 MeV, the correction factor for 
the isotopic contamination is 0.8 to 0.9. 
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FIG. 5: The experimental activation yield from ELBE 
(squares) for the ^^*Sm(7, n)^''^Sm reaction normalized to the 
photon fluence is compared to the yield calculated using cross 
sections measured by Carlos et al. (open circles). The data 
from Carlos et al. [23| have been scaled by factor of 0.8, fol- 
lowing the recommendations by Berman et al. [29| (see text). 
The downward arrow denotes the neutron emission threshold. 

For calculating the activation yield, the cross section 
data from Carlos et al. were renormalized by us to cor- 
rect for two effects: 

1. The target used by Carlos et al. had a contam- 
ination of 11.4% of other Sm- isotopes. We have 
subtracted the possible contamination using earlier 



2. As shown in a precision experiment by Berman 
et al. [29], the cross sections measured at Saclay 
for nuclei with A < 150 are larger by approxi- 
mately 20% in all cases except i^ip,. -f^p which a 
comparison was made. For this particular case, an 
additional correction factor of 0.8 was applied. 

The corrected cross sections were then used to calculate 
the activation yield which compares to our experimental 
yield as shown in Fig. [5l 

In addition, the activation yield Ym for the 
^'*'*Sm(7,n)-'^'*^™Sm reaction has been determined using 
the method described in Sec. IIII Al With known Yg and 
Y,n, the isomeric yield ratio 77 = y +y ^'"^ determined 
(see Fig.©. 

In an activation experiment by Mazur et al. [30| us- 
ing bremsstrahlung beams from the M-30 microtron facil- 
ity at Uzhgorod (Ukraine) , the isomeric-state excitations 
have been investigated for several A^ = 82 closed-shell nu- 
clei in the energy range 8-18 MeV. The (7, n)™ cross sec- 
tions were deduced using the Penfold-Leiss method f31| , 
see Fig. 1 of Ref. [30| • Mazur et al. have also measured 
the isomeric yield ratio 77 = y ™ 32 1 as a function of 

energy for the ^^''Sm(7,n)^'''^^'™Sm reaction. The exper- 
imental isomeric yield ratio from the ELBE is in good 
agreement with the data from Mazur et al, see Fig. [51 
At low energies close to the reaction threshold, the ac- 
tivation yield is influenced by the different experimental 
conditions (e.g, use of a hardener, thin/thick radiators 
etc). 



B. Model calculations in the Hauser-Feshbach 
formalism 

This section is dedicated to the comparison of ex- 
perimental results to the predictions of the two ad- 
vanced model codes TALYS and NON-SMOKER which 
are based on the Hauser-Feshbach formalism. The NON- 
SMOKER code uses the neutron optical-model potential 
by Jeukenne et al. [33| with a low-energy modification 
by Lejeune [3J|. The 7-ray strength function is taken 



from Thielemann and Arnould [35|- The low-energy 
modification of the GDR Lorentzian is by McCuUagh 
et al. 36!]. The nuclear level density implemented in the 
NON-SMOKER code is based on a global parametriza- 
tion by Rauscher et al. [37] within the back-shifted Fermi- 
gas formalism. 

The TALYS calculations presented in this paper were 
performed with the current version of the code TALYS- 
1.0. The default option of this code uses the neutron 
optical-model potential parameterizations by Koning and 
Delaroche [33] ■ The El photon strength function is from 
the compilation by Kopecky and Uhl [33|- The nuclear 
level density model is based on an approach using the 
Fermi-gas model |4q |. 

The measured activation yield (Yg) is compared to the 
simulated yield from Hauser-Feshbach models as shown 
in Fig. [71 It is observed that the experimental yield 
roughly agrees to the simulated yield using cross sections 
predicted by the default inputs to TALYS (agrees within 
20%) and NON-SMOKER (agrees within a factor of 2). 

The statistical model calculations are sensitive to the 
basic input ingredients like optical potentials, strength 
functions and level densities. In the case of (7,n) reac- 
tions, the influence of the neutron optical potential is 
weak, whereas the photon strength function is a crucial 
ingredient of the model calculations. We modified the de- 
formation dependent parameters of the El strength func- 
tion used in TALYS according to a new parametrization 
explained below (see also Ref. [4]|). 

The strength function fi{E^) is related to the average 
photoabsorption cross section 'a^{E~f) over a large num- 
ber of levels with same spin and parity by the equation 



2Jo + 1 '^{E^) 



2J^ + 1 {nhcyE^ 



MEy) 



(6) 



Jo and Jx denote the spins of ground and excited states 
respectively [43 . Eq. [6| can be applied to a wide energy 
range around the GDR, including the low energy tail of 
it. In even nuclei, the spin Jx is equal to the multipolarity 
of the GDR, i.e., A = 1. 

In nearly all types of heavy nuclei, the GDR can be 
treated as splitted into two or three components, whose 
energies are well predicted by the finite range droplet 
model (FRDM) ^. Making use of the fact that the vi- 
brational frequency Ek/h along a given axis k is inversely 
proportional to the corresponding semi-axis length Rk, 
and the splitting [4J] is due to the three different axes 
of the ellipsoid, the nuclear shape with its quadrupole 



deformation parameter /3 and triaxiality parameter 7 is 
parameterized as: 

Eq 



Eh — 



Eq ■ i?0 

Rk 



exp 



• (3 ■ cos(7 



pTT) 



(7) 



Here, Rq denotes the nuclear radius and Eq is the GDR 
centroid energy. For a spherical nucleus with mass A, Rq 
is given by Rq = 1.16^^'''^ fm and £^0 is calculated with 
an effective nucleon mass 771*= 874 MeV/c^. 

In the proposed phenomenological ansatz we parame- 
terize the strength function by a sum of three Lorentzians 
tailing down from the GDR. The centroid energy of the 
GDR is directly related to parameters of the finite range 
droplet model with one extra parameter which is equal 
for all heavy nuclei, the reduced nucleon mass m* 45]. 
As was shown theoretically [46] a description of the El 
photoabsorption cross section by a Lorentzian is appro- 
priate although the total width F of a GDR in a heavy 
nucleus is dominated by spreading and not by escape, i.e. 
direct decay. 

The average absorption cross section in the GDR is 
given by 



(7-y[E-y 



1.29-Z -N 



E^^Tk 



A 



t^El-E^^Y + EiTl 



(8) 



■^7 fc 



Here we assume that the GDR width Ffe to be used in the 
sum of the three Lorentzians depends on the resonance 
energy E^ only, in contrast to earlier descriptions |39ll47l |. 
The photon energy E^ is given in MeV and a^{E^) in 
fm^. 

In this description, the Thomas-Reiche-Kuhn sum rule 
as determined from general quantum mechanical argu- 
ments [4a] is included for obtaining the average photoab- 
sorption cross section on an absolute scale. We use re- 



sults from hydrodynamical considerations [44l | by adapt- 
ing surface dissipation to the Goldhaber- Teller model of 
the GDR. This results in a power law dependence of the 
resonance width on the respective resonance energy E)^. 
In general, for all stable nuclei with A > 80, 



TkiEk) = 0.05 -{Ekr 



(9) 



For the exponent, the value 6 — 1.6 was derived from the 
one-body dissipation model |4j| and the proportionality 
parameter stems from a fit to many nuclei including axi- 
ally as well as triaxially deformed nuclei where the GDR 
is split into three components with well defined distances 
from the centroid and relative intensities. We generalize 
this by using the proposed [4J power law dependence for 
different nuclei and thus reduce the number of parameters 
for the description of the GDR's by relating the spread- 
ing width in all nuclei with mass A > 80 to the respective 
resonance energy of their GDR. It is stressed here again 
that the Lorentzian photon strength function f{Er) is 
parameterized in each nucleus by a constant spreading 
width T{Ek) not depending on the photon energy. 



The above parametrization has been discussed in detail 
with examples from experimental studies of 8 nuclides be- 
tween A = 80 and 238 in a recent publication [4l[. The 
new phenomenological description based on the ground 
state deformation parameters describes well the electric 
dipole excitations for nuclei with A > 80 from E^ ^ 4: 
MeV up to several MeV above the GDR. In addition, the 
dipole strength function for the case of ^^^Au was deter- 
mined using the proposed parametrization and compared 
with experimental results from ELBE |15l |. 

The yield for the ^'*'*Sm(7, n) reaction calculated using 
the TALYS code with modified El photon strength func- 
tions derived on the basis of the above discussion is shown 
in Fig. [71 As in the quasi-spherical ^^''Sm the three GDR 
components coincide, the splitting is negligible and the 
following resonance parameters were used: Eq = 14.84 
MeV, ctq = 356.6 mb and Tq = 3.74 MeV. In Fig.[8l the 
isomeric yield ratio r/ measured at ELBE is compared 
to the isomeric yield ratio calculated using cross sections 
predicted by the TALYS code with default and modified 
inputs. The NON-SMOKER model does not provide the 
cross sections arising from the decay of the metastable 
state in ^^'^Sm separately and hence it was not possible 
to calculate a yield ratio. 
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FIG. 7: Experimental activation yield relative to the photon 
fluence for the ^'*'*Sni(7, n) reaction compared to theoretical 
model calculations. The experimental data are denoted by 
squares with a downward arrow denoting the neutron emission 
threshold. The dashed and dotted lines denote yield calcula- 
tions using cross sections from TALYS and NON-SMOKER 
codes respectively. The solid line represents a TALYS calcu- 
lation with modified inputs, see text. 



C. Activation yield for ^''*Sm(7,p) and ^'''*Sni(7, a) 
reactions 

The activation yields for the "'^'*^Sm(7,p) and 
^^'*Sm(7, a) reactions compared to the simulated yields 
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FIG. 8: The isomeric yield ratio 77 



Y^^Y„ 



measured at 



ELBE for the ^^^Sm(7,n) reaction (diamonds) is compared to 
the isomeric yield ratio calculated using cross sections pre- 
dicted by the TALYS code (dashed line). The solid line rep- 
resents the yield ratio calculated using modified inputs to the 
photon strength functions, see text. 
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FIG. 9: Experimental activation yield relative to the pho- 
ton fluence for the (7,p) and (7,0) reactions in ^^^Sm. The 
yields for ^**Sm (7,p) (squares) and ^^*Sm(7,Q:) (diamonds) 
are shown. The dashed and dotted lines denote yield calcula- 
tions using cross sections from TALYS and NON-SMOKER 
codes with default inputs whereas the solid line represents a 
TALYS calculation with modified inputs, see text. 



using cross sections predicted by the TALYS and NON- 
SMOKER calculations are shown in Fig. [9l The predic- 
tions given by the TALYS code with modified inputs to 
the El photon strength function is also shown. 

The uncertainties in the experimental yield shown in 
Fig. [3] are mainly the statistical uncertainties. The sys- 
tematic uncertainties arise from the photon emission 
probabilities for the decays stemming from the (7,^) 
and (7, a) (8%) reactions (see Table H]) and from 



the full-energy peak efficiency calculations (2%). 
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FIG. 10: The photoneutron cross sections a [cr(7,n)+a(7,np)] for Sm-isotopes as described_^ the photon strength parametriza- 
tion (solid lines) from Junghans ei al. [4ll ]. The original cross sections from Carlos et al. [27| have been scaled by a correction 
factor of 0.8, see text. The dashed lines denote yield calculations using cross sections from the NON-SMOKER code with 
default nuclear physics inputs. 



V. CONCLUSIONS 

Activation experiments of the kind described - espe- 
cially when performed with a transport system allowing 
the study of short half- lives deliver exact photon strength 
data. Thus they allow a verification of data obtained 
previously by direct absorption experiments or by detect- 
ing the neutrons from the process. In Fig. [TUl the pho- 
toneutron cross sections of the Sm-isotopes as described 
by the photon strength parametrization from Junghans 
et al. [4lj] is shown. 

From the figure, it is obvious that the parametrization 
describes all the isotopes well except ^'^''Sm in which the 
contribution from ^''^Sm(7,7') and ^^'*Sm(7,p) is of sig- 
nificant magnitude at low energies close-around S„ (see 
Fig. [3]). For all cross section data from Carlos et al, 
a correction factor of 0.8 based on Berman et al. [2^ 
experiments was applied. For the ^^^Sm cross sections, 
the additional correction factor arising from isotopic con- 
taminations was also necessary (see Sec. lIVAp . The 
yield calculations using cross sections from the NON- 
SMOKER code with default nuclear physics inputs are 
also shown. It is clear the calculations incorporating the 
new parametrization provides a better description to the 
experimental data than default NON-SMOKER results. 

The agreement between the parametrization and ex- 
perimental data is surprisingly good- taking into account 
that all parameters (except /3) were obtained from global 
fits to all nuclei with A > 70 (as long as respective data 



are available for them). The only non-global parameter, 
the deformation /3 was directly derived from the electric 
quadrupole transition strength to the first 2+ state p9| . 
The approximation for the triaxiality 7 is from the sys- 
tematics presented by [50] . The values of /3 and 7 used 
for each isotope are shown in the respective figure. The 
Lorentzian resonance integral corresponds to the TRK 
sum rule and the resonance energies are calculated using 
the symmetry energy J=32.7 MeV and surface stiffness 
Q=29.2 MeV [H. 

To summarize, this paper presents the photodisinte- 
gration of ^^"^Sm to four exit channels. The ^''*Sm(7,n) 
reaction yield has been compared to the yield calculated 
using cross sections from previous photoneutron experi- 
ments and a comparison of the two data sets leads to a 
conclusion on the inaccuracies in previous data. It has 
been verified that the recommendations based on Berman 
et al. [29] experiments are necessary to correct the cross 
section measurements at Saclay. 

The activation yield for all the photodisintegration 
reactions has been compared with TALYS and NON- 
SMOKER models. In general, the experimental activa- 
tion yields agree within a factor of 2 to the simulated 
yield using statistical model predictions. In view of the 
fact, that the codes NON-SMOKER and TALYS are 
based on the same theory-the Hauser-Feshbach model, 
eventual differences appearing in the model predictions 
are strongly related to their nuclear physics inputs. The 
effect of using different photon strength functions in the 
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TALYS model calculations has been discussed here. For 
the ^^^Sm(7,Q;)reaction, the statistical model description 
faces the additional problem of the possible emission of 
low energy photons prior to the alpha emission. In a 
forthcoming paper the respective data for ^'*^Sm will be 
presented together with results for a different mass re- 
gion, and the proper treatment of that process by the 
two codes will be examined. 

The new proposed phenomenological parametrization 
describes very well the photoneutron cross section data 
for Sm isotopes. A similar investigation has been per- 
formed for the stable even isotopes of Mo and same effects 
has been observed J51|] ■ It is thus worthwhile to point out 
again that this new fit only needs a very limited number 
of fit parameters to describe the electric dipole strength 



in all nuclei with ^ > 70, for which nuclear photoeffect 
data exist. 



VI. ACKNOWLEDGEMENTS 

We thank P. Michel and the ELBE team for provid- 
ing a stable beam during activation experiments. We are 
indebted to J. King for the MCNP simulations of the 
bremsstrahlung spectra at the photoactivation site. Spe- 
cial thanks are due to M. Kohler and D. Degering, for 
their help with the measurements at the underground 
laboratory Felsenkeller, Dresden. The technical assis- 
tance of A. Hartmann is gratefully acknowledged. 



[1] E. M. Burbidge, G. R. Burbidge, W. A. Fowler et al, 

Rev. Mod. Phys. 29, 547 (1957). 
[2] D. L. Lambert, Astron. Astrophys. Rev. 3, 201 (1992). 
[3] M. Arnould and S. Goriely, Phys. Rep. 384, 1 (2003). 
[4] S. S. Dietrich and B. L. Berman, At. Data Nucl. Data 

Tables 38, 199 (1988). 
[5] J. Audouze and D. N. Schramm, Nature 237, 447 (1972). 
[6] R. Andreasen and M. Sharma, Science 314, 806 (2006). 
[7] E. Somorjai, Z. Fiilop, A. Z. Kiss, et al., Astron. Astro- 
phys. 333, 1112 (1998). 
[8] W. Rapp, J. Gorres, M. Wiescher, et al., Astrophys. J. 

653, 474 (2006). 
[9] G. Rusev, R. Schwengner, F. Donau et al., Phys. Rev. C 

73, 044308 (2006). 
[10] R. Schwengner, G. Rusev, N. Benouaret et al., Phys. Rev. 

C 76, 034321 (2007). 
[11] R. Schwengner, G. Rusev, N. Tsoneva et al., Phys. Rev. 

C 78, 064314 (2008). 
[12] G. Rusev, R. Schwengner, F. Donau et al., Phys. Rev. C 

77, 064321 (2008). 

[13] G. Rusev, R. Schwengner, R. Beyer et al., Phys. Rev. C 

79, 061302 (2009). 
[14] N. Benouaret, R. Schwengner, G. Rusev et al., Phys. Rev. 

C 79, 014303 (2009). 
[15] C. Nair, M. Erhard, A. R. Junghans et al, Phys. Rev. C 

78, 055802 (2008). 

[16] R. Schwengner, R. Beyer, F. Donau et al., Nucl. lustrum. 

Methods A 555, 211 (2005). 
[17] A. Wagner, R. Beyer, M. Erhard et al, J. Phys. G: Nucl. 

Part. 31, S1969 (2005). 
[18] MCNP - Monte Carlo N-Particle Transport Code, 

http:/7mcnp-green. lanl.gov/; Version 4C2 obtainable 

from http://www.nea.fr/abs/html/ccc-0701.html 
[19] S. M. Seltzer and M. J. Berger, At. Data Nucl. Data 

Tables 35, 345 (1986). 
[20] C. Nair, A. R. Junghans, M. Erhard et al, J. Phys. G; 

Nucl. Part. 35, 014036 (2008). 
[21] M. Kohler, D. Degering, M.Laubenstein et al., Appl. Rad. 

Isot. 67, 736 (2009). 
[22] Evaluated Nuclear Structure Data Files (ENDSF). The 

database service of National Nuclear Data Center 

(NNDC) . 'http://www.nndc.bnl.gov/ensdf/ 
[23] GEANT3: CERN program Library Long Writeup Q121, 



CERN, Geneva (CH), 1994. 

[24] A.J. Koning, S. Hilaire and M.C. Duijvestijn, Pro- 
ceedings of the International Conference on Nuclear 
Data for Science and Technology, April 22-27, 2007, 
Nice, France, editors O.Bersillon, F.Gunsing, E.Bauge, 
R.Jacqmin, and S.Leray, EDP Sciences, p. 211-214 
(2008); http://www.talys.eu 

[25] T. Rauscher and F. -K. Thielemann, At. Data Nucl. Data 
Tables 88, 1 (2004). 

[26] T. Vidmar, Nucl. lustrum. Methods A 550, 603 (2005). 

[27] P. Carlos, H. Beil, R. Bergere, et al, Nucl. Phys. A 225, 
171 (1974). 

[28] R. Bergere, H. Beil, P. Carlos, et al, Nucl. Phys. A 133, 
417 (1969). 

[29] B. L. Berman, R. E. Pywell, S. S. Dietrich et al, Phys. 
Rev. C 36, 1286 (1987). 

[30] V. M. Mazur, V. A. Zheltonozhskii, and Z. M. Bigan, 
Phys. of Atomic Nuclei 58, 898 (1995). 

[31] A. S. Penfold and J. E. Leiss, Phys. Rev. 114, 1332 
(1959). 

[32] R. Vanska and R. Rieppo, Nucl. Instrum. Methods 179, 
525 (1981). 

[33] J. P. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rev. 
C 15, 10 (1977). J. P. Jeukenne, A. Lejeune, and C. Ma- 
haux, Phys. Rev. C 16, 80 (1977). J. P. Jeukenne, A. 
Lejeune, and C. Mahaux, Phys. Rep. 25, 83 (1976). 

[34] A. Lejeune, Phys. Rev. C 21, 1107 (1980). 

[35] F. -K. Thielemann, and M. Arnould, Proc. Int. Conf. 
Nucl. Data for Science and Technology, edited by K. 
Bockhoff, D. Reidel Publ. Co, Brussels, 762 (1983). 

[36] C. M. McCuUagh, M. L. Stelts, and R. E. Chrien, Phys. 
Rev. C 23, 1394 (1981). 

[37] T. Rauscher, F.-K. Thielemann, and K. L. Kratz, Phys. 
Rev. C 56, 1613 (1997). 

[38] A. J. Koning and J. P. Delaroche, Nucl. Phys. A 713, 
231 (2003). 

[39] J. Kopecky and M. Uhl, Phys. Rev. C 41, 1941 (1990). 

[40] A. J. Koning, S. Hilaire and S. Goriely, Nucl. Phys. A 
810, 13 (2008). 

[41] A. R. Junghans, G. Rusev, R. Schwengner et al, Phys. 
Lett. B 670, 200 (2008). 

[42] G. A. Bartholomew, E. D. Earle, A. J. Ferguson et al. 
Adv. Nucl. Phys. 7, 229 (Plenum Press, New York, 1973). 



11 



[43] W. D. Myers, W. J. Swiatecki, T. Kodama et at, Phys. 

Rev. C 15, 2032 (1977). 
[44] B. Bush and Y. Alhassid, Nucl. Phys. A 531, 27 (1991). [49] 

[45] P. MoUer, J. R. Nix, W. D. Myers et ai, At. Data Nucl. 

Data Tables 59 185 (1995). [50] 

[46] C. B. Dover, R. H. Lemmer and F. J. W. Hahne, Ann. 

Phys. 70 458 (1972). [51] 

[47] L. Zanini, F. Corvi, H. Postma et ai, Phys. Rev. C 68, 

014320 (2003). 
[48] J. M. Eisenberg and W. Greiner, Nuclear Theory Vol. 



2: Excitation mechanisms of the nucleus, Third edition. 

North Holland Physics Publ., Amsterdam (1988). 

S. Raman, C. W. Nestor, Jr., and P. Tikkanen, At. Data 

Nucl. Data Tables 78, 1 (2001). 

W. Andrejtschoff and P. Petkov, Phys. Rev. C 48, 2531 

(1993). 

M. Erhard, A. R. Junghans, C. Nair et ai, Phys. Rev. C 

81, 034319 (2010). 



